Do Vortices Entangle? 
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We propose an experiment for directly constructing and locally probing topologically entangled 
states of superconducting vortices which can be performed with present-day technology. Calculations 
using an elastic string vortex model indicate that as the pitch (the winding angle divided by the 
vertical distance) increases, the vortices approach each other. At values of the pitch higher than a 
maximum value the entangled state becomes unstable to collapse via a singularity of the model. We 
provide predicted experimental signatures for both vortex entanglement and vortex cutting. The 
local probe we propose can also be used to explore a wide range of other quantities. 



PACS: 74.25.Qt, 74.25.Sv 

The high superconducting transition temperatures 
of compounds such as YBa 2 Cu30( 7 _5) (YBCO) and 
BiSr^CaaCuOg (BSCCO) lead to a very rich set of behav- 
iors of the magnetic vortices which form inside the mate- 
rial in the presence of an applied magnetic field. Thermal 
fluctuations are significant over a wide range of the (H,T) 
phase diagram [1] , causing the lattice of stiff vortex lines 
to melt well below T c . The nature of this molten state 
has remained a subject of intense debate. 

Since there can be significant thermally-induced wig- 
gling along the length of the vortex in the liquid state, 
Nelson proposed that neighboring vortex lines may be- 
come entangled with each other, in analogy with a super- 
fluidity transition in a boson system [2]. This entangle- 
ment could produce a dramatic increase in the viscosity 
of the vortex liquid [3-5] , similar to that which occurs for 
entangled polymers. As a result, vortex pinning by ran- 
dom disorder in the sample would be enhanced, so that 
the resistivity would drop in the entangled state. 

In order for the vortices to entangle, it is crucial that 
neighboring lines not cut through each other easily [6] 
and reconnect into a disentangled state. Estimates of 
the cutting barrier vary widely [3,5,7-10], ranging from 
50ksT to of order ksT, leaving the question of whether 
vortices can easily cut in the liquid state unresolved. Nu- 
merical simulations performed in limits ranging from the 
low-field London regime to the high-field lowest Landau 
level regime have proven similarly ambiguous, with some 
simulations interpreted as providing evidence for entan- 
glement [11] and others interpreted as showing that the 
lines cut and do not entangle [9,12,13]. The simulations 
are limited both by the models chosen and by the system 
sizes that can be simulated. Models based on the boson 
analogy lack long-range interactions along the vortex line, 
which can stiffen the vortices and might reduce entangle- 
ment. [7,14,15]. In the frustrated 3D XY model, there 
are multiple ways to define a path followed by a vortex 
line, some of which are consistent with entanglement, and 
others which are not [12]. 

Since theoretical and numerical evidence for entangle- 
ment have proven inconclusive up to this point, it is nat- 



ural to turn to experiments to resolve the issue. Unfortu- 
nately experimental evidence for or against entanglement 
[16-19] has also proven ambiguous, both because the ex- 
perimental measures are indirect, and because there is 
considerable variation in the definition of entanglement 
[20], such as equating entanglement with the onset of 
plasticity [16]. Thus, despite more than a decade of theo- 
retical, numerical, and experimental studies, the question 
of whether vortices in high-temperature superconductors 
can form an entangled state has not yet been convinc- 
ingly answered. 

Here we propose a direct experimental test of vortex 
entanglement by means of a local atomic-force micro- 
scope (AFM) probe which can unambiguously determine 
whether it is possible for two vortices to wind around 
each other without cutting. We consider conditions that 
are as favorable as possible for entanglement: low vortex 
density, weak background pinning, and low temperatures, 
which should increase the barrier to flux cutting. We pro- 
pose using a magnetic AFM tip to produce the simplest 
possible entangled state, two vortices wound around each 
other, and show theoretically that the forces involved in 
creating the entangled state can be measured with cur- 
rent technology. Should the vortices cut rather than en- 
tangling, a distinct force signature will be observed. If en- 
tanglement does not occur in this limit, it seems unlikely 
to occur in the higher density, higher temperature condi- 
tions near the melting transition. In addition to testing 
entanglement, the local AFM probe which we propose 
can be used to measure a wide range of other quantities, 
such as vortex line tension, local pinning force, and local 
shear forces. 

We use as our definition of an entangled vortex state 
the picture originally proposed by Nelson, in which vor- 
tices behave like elastic lines that cannot easily cut 
through each other [2]. We consider a layered super- 
conductor containing two closely spaced magnetic dots 
at the bottom of the sample [21] (dark circles in Fig. 1) 
A mobile magnetic dot is introduced to the top of the 
sample in the form of a magnetic AFM tip (open circle 
in Fig. 1). An optional additional magnetic dot 
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FIG. 1. (a) Schematic of the starting state. Black dots 
represent the fixed pins; the open circle represents the end 
of the vortex being moved by the AFM. (b) A configuration 
of wound vortices, found numerically by minimizing the free 
energy within an elastic string model. The lengths are in units 
of A and the system is isotropic. 

can be added at the top of the sample to pin one of the 
vortices. We assume that there is a sufficiently small ex- 
ternally applied magnetic field that there are only two 
vortex lines within the region of interest, an area a few 
A on a side, where A is the London penetration depth 
of the material. The magnetic dots attract the vortices, 
so that the top and bottom positions of the vortices are 
fixed at the dot locations. As the AFM is moved along 
the surface of the sample, it drags the top of one vortex. 
The force required to drag the vortex line can be di- 
rectly measured. When the AFM tip moves in a circular 
path around the fixed upper dot, the two vortices wind 
together, producing an entangled state. In Fig. 1(b), we 
illustrate one possible entangled vortex configuration, ob- 
tained numerically as described below. As the winding 
angle increases, the force required to drag the vortex fur- 
ther around increases according to a form derived below. 
If the vortices cut, this force will abruptly drop. Thus, 
using such an experiment, it is possible to directly probe 
whether vortex entanglement can occur. 

To estimate the force required to entangle a pair of 
linclike vortices, we represent the vortex lines as one- 
dimensional elastic strings that cannot cross. Such a 
model may be applicable in the London limit at low fields 
B < 0.2B c2 , although the lack of long-range interactions 
along the z axis is a notable limitation [7,14,15]. The 
vortex configuration is determined by a balance of the 
interaction energy, which drives the strings apart, and 
the elastic energy, which pulls the strings together to re- 
duce their length as they are wound. Thus, the greater 
the pitch (defined as the rate of change of winding angle 
per vertical distance), the closer the vortices approach, 
as they seek the preferred spacing that minimizes their 
free energy. Even if they are held radially away from the 
preferred spacing, as is shown at the top of Fig. 1(b), 
they approach the preferred spacing within the middle 



of the sample, leading to the roughly constant spacing 
in the middle of the sample as shown. This applies up 
to a certain pitch; beyond that pitch, we will show that 
the vortices are unstable to a collapse. A similar insta- 
bility in a related model was discussed in Ref. [8]. The 
existence within the model of a collapse of the entangled 
state at a singularity when the vortices approach each 
other too closely raises the question of whether the en- 
tangled state can exist in other experimental situations 
where the vortices are driven closer than this distance. 

Elastic String Model — We consider two vortices, at 
positions ri(z), ^(z), with r = (x,y), giving the position 
in the plane as a function of the vertical distance z. As 
a starting point, we consider the free energy from the 
elastic string model [3]: 

L 

F = J dzY,\h{d z r l ) 2 + 2Y,eM\r l {z)~r 3 {z)\/X), 

i i<3 

(1) 

with e = ($o/4ttA) 2 , l x m y/M ± /M z e . We fix the 
boundary conditions on the vortex position at the top 
(z = L) and the bottom (z — 0) of the sample, and 
then minimize the free energy (1) to find the positions 
of the vortex in between. To specify these positions, 
we need eight real numbers, two for each vortex at the 
top and another two for each vortex at the bottom. To 
simplify, let us fix two of these coordinates by setting 
fl(0) + r*2(0) = 0; thus, the origin of the coordinate sys- 
tem is set at the midpoint of the two vortices at the 
bottom of the sample. 

We note that if a given pair of functions rx^(z) min- 
imize the free energy for given boundary conditions 
fi^(L), then the functions fi. 2(2) + vz also minimize the 
free energy for a different set of boundary conditions: 
ri,2{L) + vL. Thus, it suffices to consider only the case 
in which also r\(L) + f^iL) = 0, as then all other bound- 
ary conditions at the top can also be obtained. In this 
case, for all z, n(z) + T2[z) — 0. Thus, introduce co- 
ordinates with r\(z) — \r[z) cos(0(z)), r(z) sin(0(z))] and 
ri{z) = —ri(z). For the boundary conditions, we fix 
r(0), r(L), 8(0), 0(L). Then, the extremization of the free 
energy yields equations of motion: 

r 2 9 z = J (2) 

h{d 2 r - J 2 /r 3 ) + (2e Q )K 1 (2r/X)/X = 0. (3) 

Exploiting an analogy of this system to a particle evolv- 
ing in time, the first equation is the familiar conservation 
of angular momentum. 

There exist solutions with constant vortex spacing, 
r(z), and pitch, d z 9(z). To find these, set d 2 r = in 
the equations of motion to obtain 
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FIG. 2. Angular force on the top end of each vortex as a 
function of pi tch d z 9. Force is measured in units of ei, pitch 
in units of ^/eo/ei/A. The dashed line denotes an unstable 
vortex configuration above a pitch of approximately 0.378, 
corresponding to a closest vortex spacing of r « 1.19A. 

This gives the preferred pitch ("IS M function of spacing, 
by a balance between elastic energy of one vortex and 
interaction energy between two. For small r, this reduces 
to d z 9 = (l/r)y / e u /?i. From this one can show that 
for small spacing in the isotropic system (eo = ei) the 
vortices cross at an angle of ir/2, as expected [8]. While 
this scenario seems reasonable, we will find later that 
there is a major caveat: Eq. (4) describes a minimum of 
the free energy for large r, but for smaller r it describes 
only an extremum and is unstable to a collapse of the 
vortices. 

To find the force that would be detected by the AFM 
tip as the vortices are wound around one another into 
an entangled state, we first consider the large r case of 
Eq. (4) when the entwined vortex configuration is sta- 
ble. We fix the total winding angle, A9 = 9(L) - 6(0), 
and the values r(0),r(L). The force that the dragged 
vortex exerts on the AFM tip can be found by taking a 
derivative of the free energy with respect to r and 9. We 
find that the angular force on the top end of each vortex 
is eird z 9(L), while the radial force on each is e\d z r(L). 
The radial force will vanish when r, d z 8 obey Eq. (4) at 
the top and bottom of the sample, giving a solution with 
constant radius and pitch d z = (9(L) — 9(0)) /L. In this 
case, we plot the angular force as a function of pitch in 
Fig. 2. 

Experimental Implications- The magnitude of the an- 
gular force is large enough to be detected experimentally. 
For example, consider two vortices in a YBCO sample. 
For this material, l\ ~ 5eo, where eo ~ 140 pN. Thus, 
in Fig. 2, a force of 0.1 l\ would correspond to approx- 
imately 70 pN, a value within the range of forces de- 
tectable with AFM. In order to convert the pitch d z 9 
plotted in Fig. 2 into the total angular displacement 
A9 imposed on the dragged vortex, we must know the 



thickness L of the sample in the z direction. Assuming 
L = 1/rni gives AO = Ld z 9 = 0.91tt for d z 9 = 0.2. The 
spacing of the pins at the bottom of the sample does not 
affect the force measured at the top of the sample unless 
it is of order the sample thickness or wider, due to the 
fact that the vortex spacing in the middle of the sample 
depends only on the pitch. 

If the radial force is non-vanishing, then r(z) is not a 
constant function: by applying a radial force at the top 
we drive the vortices away from the preferred spacing. If 
r(L) is too large for the pitch, then r(z) will decrease in 
the middle of the sample, while if r(L) is too small, r(z) 
will increase in the middle of the sample. Solutions to the 
equations of motion can be found numerically [22], as in 
Fig. 1. For large enough L, one finds that in the middle 
of the sample the r, d z 9 are given with good accuracy by 
Eq. (4), for large enough r when the equation describes 
a minimum. 

We next consider under what conditions the twisted 
vortex configuration becomes unstable. This occurs when 
Eq. (4) no longer describes a minimum. We find by differ- 
entiating the free energy twice that for stability at fixed 
J, we need 

372 o 

h— - e — (K (2r/\) + K 2 (2r/X)) < 0. (5) 

When Eq. (5) vanishes, at r w 1.19A, the system is 
marginal, and for smaller r, the system is unstable to 
perturbations. In this case, the repulsion between the 
vortices is unable to overcome the elastic energy, and the 
two vortices will be driven together, meeting at a singu- 
larity for some z with r(z) = 0. Only a finite energy cost 
is paid to have the two vortices meet at a point, while 
an arbitrary amount of winding can be accomplished at 
that point with no free energy cost. Within an elastic 
string model, it can be shown that for any vortex-vortex 
interaction energy which diverges less strongly than 1 /r 2 
there is an instability of this nature for sufficiently small 
r. At the marginal r « 1.19A, the angular force on each 
vortex is equal to (0.45...)v / eieo- Thus, this is the largest 
angular force possible in this experimental setup, while 
the minimum stable vortex spacing is 2r w 2.39A. 

In more highly anisotropic materials such as BSCCO, 
the elastic string model for the vortex lines is expected to 
break down due to the possibility of decoupling of vortex 
pancakes on adjacent layers. Consider a magnetic AFM 
tip bound to a pancake at the top of a vortex line in 
BSCCO. If this pancake decouples from the rest of the 
vortex line below, then the angular force as a function 
of angular displacement will not follow the form shown 
in Fig 2. but will instead be periodic in the winding 
angle, due to the interaction with the remaining stack 
of pancake vortices left behind. This is illustrated by a 
simulation of a pancake vortex system [23] shown in Fig. 
3. The force from a single pancake is extremely 



3 




Winding Angle 



FIG. 3. Angular force on a decoupled pancake vortex in 
units of ?i as a function of total winding angle A9, obtained 
from a simulation of pancake vortices. 

small, of order 10~ 20 N; thus it is natural to ask how to ex- 
perimentally determine whether one pancake is attached 
to the AFM tip, or whether the entire vortex stack has 
detached from the AFM tip so that the tip is moving 
freely over the sample. The presence of a vortex under 
the AFM tip can be detected by means of a local density 
of states measurement performed by temporarily chang- 
ing the mode of operation of the tip to a tunneling probe. 

Discussion — The local experimental probe that we 
propose can also be used to explore numerous other prop- 
erties of the vortex system besides vortex entanglement. 
For example, in a geometry containing only one magnetic 
pin and one vortex line, the AFM tip can be used to mea- 
sure the vortex line tension directly. If the line tension 
is known, the tip could be used to tear a vortex away 
from a (weaker) individual pin, such as a grain boundary, 
and the pinning force could be measured. By applying 
a transport current to the sample, the Lorentz force can 
be determined directly. Local rheology measurements are 
also possible in the vortex lattice state; for example, the 
local elastic constants can be probed by moving a single 
vortex back and forth around its lattice equilibrium po- 
sition. The temperature dependence of both the elastic 
constants and the pinning energy could also be probed. 

Conclusion — We have proposed an experimental setup 
for constructing and probing entangled states of super- 
conducting vortices, and shown that the forces associ- 
ated with vortex entanglement are experimentally mea- 
surable. This kind of experimental setup can be gener- 
alized to other types of vortices, as in fluid turbulence. 
Within the elastic string model, we find that the entan- 
gled state is only stable up to a maximum pitch or min- 
imum vortex spacing. The instabilities we have found 
raise the question of whether the entangled state can ex- 
ist with a high density of vortices. To answer this ques- 



tion and to compare pancake and elastic string models, 
an experimental test of our proposal is desirable. 
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